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Abstract. - Realistic charged macromolecules are characterized by discrete (rather than 
homogeneous) charge distributions. We investigate the effects of surface charge discretization 
on the counterion distribution at the level of mean-field theory using a two-state model. Both 
planar and cylindrical geometries are considered; for the latter case, we compare our results 
to numerical solutions of the full Poisson-Boltzmann equation. We find that the discretization 
of the surface charge can cause enhanced localization of the counterions near the surface; for 
charged cylinders, counterion condensation can exceed Oosawa-Manning condensation. 



Introduction. - The interaction between charged macromolecules and counterions is a 
crucial component of the physics of charged systems [1] . Although realistic macromolecules 
are often composed of discrete charges, theoretical models typically assume a homogeneous 
surface charge distribution. In the case of cylindrical polyelectrolytes (PEs), this assumption 
leads to Oosawa-Manning (OM) condensation [2,3], in which counterions become closely as- 
sociated with the macromolecule, effectively lowering its overall charge. A similar effect also 
occurs for spherical macromolecules [4]. Theoretically, the effects of surface charge inhomo- 
geneities on the counterion distribution has been studied for planar [5,6], cylindrical [7-10], 
and spherical [11,12] geometries. In particular, both Moreira et al. [5] and Lukatsky et al. [6] 
have shown analytically that the heterogeneity of a planar charge distribution leads to an 
enhanced localization of the counterions near the macromolecular surface. Thus, theoretical 
models that assume homogeneous surface charge distributions - including OM condensation 
theory - must be modified for realistic charge distributions. 

In this letter, we explore the effects of surface charge discretization on the counterion 
distribution in the absence of added salt. For small charge modulation, it is appropriate to 
use a perturbative expansion of the Poisson-Boltzmann (PB) equation to describe the mean- 
field behavior of the counterion distribution. It has already been shown that such an expansion 
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leads to enhanced counterion localization for planes [6] , and it is straightforward to show that 
the same conclusion holds for surfaces of arbitrary geometry [13]. For a discrete surface 
charge distribution, however, pair association between discrete macromolecular charges and 
the discrete counterions [14, 15] can become important; that is, a large number of counterions 
can become localized near the macromolecular surface. In PB theory, the counterions are 
treated as a smooth distribution function, rather than as discrete objects. In calculating 
the electrostatic energy, such a description clearly fails when many counterions are spatially 
localized, because it integrates over these localized regions. As a result, the (coarse-grained) 
self-energy of the localized counterions is incorrectly included in the electrostatic energy. Thus, 
PB theory cannot sufficiently describe the counterion distribution of a discretely charged 
surface in this regime. 

Motivated by these considerations, we introduce a two-state model (TSM) which explic- 
itly allows for the possibility of surface charge-counterion pair association. Counterions in 
the bound state are confined to the surface by the electrostatic attraction to a single discrete 
surface charge; free counterions, on the other hand, "see" the macromolecule and the bound 
counterions as a homogeneous charge distribution. We ignore any correlations in the ther- 
mal fluctuations of the counterions, but explicitly take into account the discreteness of the 
bound counterions. Our main result is that, for both planar and cylindrical geometries, the 
TSM predicts that surface charge discretization leads to enhanced counterion localization. In 
particular, it is possible to have counterion condensation on PEs that exceeds the amount 
predicted by OM condensation theory. We also find that, for cylindrical geometries, the TSM 
yields excellent agreement with numerical solutions to the PB equation when the number of 
bound counterions is small, but that the PB equation dramatically under-predicts the extent 
of counterion localization when this number is large. 

Two-state model. - For a macromolecule with an arbitrary surface charge distribution, 
the mean-field counterion density profile is determined by the competition of the electrostatic 
attraction to the macromolecule with both the counterion entropy and the interaction between 
counterions. For a discrete surface charge distribution, the localized counterions will expe- 
rience a strong electrostatic attraction to a discrete surface charge. In order to capture the 
effects of surface charge discretization on the counterion distribution, we allow the counterions 
to exist in one of two possible states: in the bound state, counterions bind to a single charge on 
the surface, sacrificing most of their entropy to gain electrostatic energy; in the free state, by 
contrast, counterions gain less electrostatic energy but have a much larger entropy. Thus, the 
two-state model is essentially a picture for counterion association with, and the consequent 
charge regulation of, a discretely charge surface. The possibility of counterion association for 
homogeneously charged objects has been considered by several other authors [16-18]. 

Consider an aqueous solution of charged macromolecules and their neutralizing counteri- 
ons. We use the cell model, in which each macromolecule and its neutralizing counterions (of 
valency Q) are enclosed in a cell of some volume V; interactions between cells are ignored. We 
also work exclusively within the primitive model, where the aqueous solvent is approximated 
as a homogeneous medium with a constant dielectric constant e (e ss 80 for water). The 
surface S s of each macromolecule has a total area A s and is decorated with a lattice (with 
lattice constant d) of N s discrete surface charges of valency Q s . Let Nf be the number of 
counterions in the free state and Nf, the number of counterions in the bound state; by charge 
neutrality, Q S N S — QNf + QNb- The bound counterions lie on a surface (of total area 
Ab), which is of the same geometry as S s but separated from it by 5, the distance of closest 
approach (set, for example, by the size of the counterions). Because of electrostatic repulsion 
and excluded volume interactions, we only allow one counterion per binding site. Defining 
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the free counterion density rif(x) = J^fJi ~ an d assuming S -C d, the electrostatic 
energy can be written as 



®n f (af) + QN b n b (x') - Q s N s n s (x') 



l B QQ s N b 



v 

, ,r ,r ff d 3 xd 3 x' . , JflQ 2 iV h 2 ff d 3 xd 3 x' 

-l B QQ s N b N s JJ — — n fc (x)n s (:r') + 2 b yy |f _ f/| n 6 (g)n t (^) 

with ?s = e 2 /eksT being the Bjerrum length = 7.1 A in water), e the elementary unit 
of charge, fcs Boltzmann's constant, T the temperature, and f3 = 1/ksT. n s (x) and n b (x) 
are the distribution functions for the surface and bound charges, respectively. Note that the 
dependence of (3TL on N s and N b has been written out explicitly, so that the normalizations 
of n s {x) and n b (x) are set to unity. 

In order to estimate the terms involving the surface charges and bound counterions in 
eq. (1), we make a series of approximations. First, the bound counterions are treated at 
the level of mean-field theory. This is equivalent to distributing the charge of the bound 
counterions equally among each binding site. In addition, the discrete nature of the surface 
charges and the bound counterions can be ignored when calculating most of the terms in 
eq. (1). In the limit { « d, charge discreteness clearly must be taken into account when 
calculating the interaction of a surface charge with the counterion charge in its binding site; 
we have therefore written this term out explicitly in the electrostatic energy [the second term 
in eq. (1)]. We can approximate the other terms in eq. (1), however, by smearing out the 
surface charges and the bound counterions on S s and T, b , respectively. Consequently, we set 
n s {x) = l/A s on E s and elsewhere; similarly, we set n b (x) = 1/A b on ~E b and elsewhere. 
Since the interactions of every surface charge with the counterion charge in its binding site are 
accounted for in the second term of eq. (1), we must be careful not to also include the coarse- 
grained approximation to this interaction. Therefore, we exclude the appropriate region R\ 
from the second integral in eq. (1), so that \x — x'\ does not run over the regions in which 
this interaction is being calculated. Finally, the last integral in eq. (1) is an approximation to 
the sum IbQ 2 N 2 /(2N s ) X^i^y Vl^ — | , where &, are the binding site locations. Because this 
sum excludes the i = j terms, the integral should exclude the region R2 in which \x — x'\ lies 
within the Wigner-Seitz cell of a single binding site. 

The free energy (3F = (37i — Sf — S b , where Sf — — f d 3 xnf(x)(hi[nf(x)vo} — 1) and 
S b = \n[N s \/N b \(N s — N b )l] — N b ln[w /w&] are the entropies of the free and bound counteri- 
ons, respectively, with v n being the volume of a counterion and v b the volume within which 
each bound counterion is confined. The two-state counterion distribution is found by min- 
imizing the grand free energy, (3G = (3F — /j,[N b + J rif(x)d 3 x} with respect to nf(x) and 
N b (the chemical potential fi acts as a Lagrange multiplier which enforces charge neutral- 
ity). Minimization of G with respect to nf(x) yields the expected Boltzmann distribution, 
rif(x) = no exp [— y(x)} , with no = c^/vq. The potential y(x) obeys the Poisson-Boltzmann 
(PB) equation, V 2 y(x) = — n 2 cxp[— y(x)], where k 2 = ^ttIbQ 2 ^. Because we have ap- 
proximated n s {x) and n b (x) as homogeneous functions in eq. (1), we can use Gauss' Law to 
show that the boundary condition for y(x) at is identical to the boundary condition for a 
"renormalizcd" charge QN b — Q S N S distributed uniformly across the surface Y^ b . Thus, the 
free counterion charge distribution in the TSM is given by the solution to the PB equation for 
the appropriate geometry, but with a renormalized and uniform surface charge distribution. 

The number of bound counterions N b is found by minimizing the grand free energy with 
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respect to Nb, 

-ln(j±-l)+J d 3 xn b (x)y(x) - l B Q 2 N b J f n b (x)n b (x') (2) 

V R 2 

+IbQQ s N 8 J J -j^—^n b (x)n s (x ) - In (n w b ) — = 0. 

Ri 

The entropy of the bound counterions gives rise to the first term in eq. (2), a term which 
clearly diverges as N b — * 0. Therefore, the solution to cq. (2) must be some value of Nb > 0; in 
other words, a discretely charged surface of arbitrary geometry will always have some bound 
counterions in the two-state model. 

In order to quantify the effects of surface charge discretization on the counterion distri- 
bution, we can compare the results of the TSM with the PB solution for the corresponding 
homogeneously charged surface. To do so, we need the PB potential y(x) for a homogeneously 
charged surface of the appropriate geometry. Analytically, the solution to the PB equation is 
known only in planar and cylindrical geometries, so we restrict ourselves to these geometries 
in what follows. 

For planar macromolecules, S s is the plane z = —S, £& is the plane z = 0, and the 
enclosing cell is the positive half-space z > 0. Consider the fraction <p of counterions enclosed 
between E& and a plane at z > 0: 4>(z;Nb) = 1 — QA/(Q S N S ) f°° dz' nf(z'; N b ), where A is 
the area of the plane. For the homogeneously charged plane, the enclosed fraction is given by 
<f>(z; N b = 0); for a discretely charged plane, the number of bound counterions is determined 
by eq. (2) . If we assume that the surface charges are in a square lattice with lattice constant 
d, then both i?i and i?2 are the region {\x — x'\ < d,\y — y'\ < d}. In the limit 5 <C d, eq. (2) 
reduces to 



o — o r Qcr 



°r(l + ^) =2^Q^ eXP ( v —) (3) 



where a = Q S N S /A is the bare charge density of the plane and oy = (Q S N S — QN b )/A is the 
renormalized charge density of the plane. The free counterion density is given by nf(z; N b ) = 
[2itIbQ 2 (z + C(cr r )) 2 ] -1 , where ((a) = 1/(2ttIbQct) is the Gouy-Chapman length [1]. Clearly, 
rif(z;N b ) is a monotonically decreasing function of N b for any z > 0. This in turn implies 
that (f>(z; N b > 0) > cf>(z; 0) for any z > 0. Thus, the two state model predicts that charge dis- 
cretization of a planar macromolcculc leads to enhanced counterion localization, in agreement 
with the results of [5,6]. 

For cylindrical macromolecules, we consider the simplest cylindrically symmetric charge 
distribution, a linear array of charges separated by a distance d (the results are qualitatively 
unchanged for a charged cylinder with a non-zero radius). In this case, S s is the line z = 0, 
Et is the surface of a cylinder of radius <5, the enclosing cell is a cylinder of radius R, and both 
i?i and i?2 are the region \z — z'\ < d. cf) is defined to be the fraction of counterions enclosed 
between T, b and a cylinder of radius S < r < R: <f>[r\ Nb) = l — 2irLQ/(Q s N s ) dr' n f (r' ; N b ) , 
where L is the length of the cylinder. Again, the enclosed fraction is given by <j>(r; N b = 0) 
for the homogeneously charged cylinder; for the discretely charged cylinder, eq. (2) becomes, 
in the limit i<d, 



1 Qs 



(7r + (l-€r))(^J =— «P( — )=<* W 
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where £ = lbQQ s /d is the bare Manning parameter, £ r = IbQ(Q s N s — QNb)/L is the renor- 
malizcd Manning parameter , £o = p/0- + p)i P = ln[-R/£], an d the integration constant j r is 
given by the implicit equation £ r = (1 + 7r)/(l + 7r cot [7 r p]), with < -f r < n/p. Eq. (4) 
applies for £ r > £o; for £ r < £o, 7r — » «7r, with < 7 r < 1. The crossover £ r = £o occurs 
when a = a = (1 + p)~ 2 [(£ - Co) -1 - Qs/Q£](<V 2c 250 ; for a < ao, (, r < £o, and vice-versa. 
Note that, since eq. (4) always gives £ r < £, this crossover can occur only when £ > £o- For 
cylindrical geometries, the free counterion density is given by [19, 20] 



for £ r > £o (again, for £ r < £o, 7r - * *7r)- Careful examination of eq. (5) reveals that, for any 
r, n/(r) is a monotonically increasing function of £ r (i.e. a monotonically decreasing function 
of Nb). In other words, 4>(r:Nb > 0) > <p(r; 0) for any S < r < R. Thus, as in the case of 
planar geometries, we conclude that charge discretization of a cylindrical charge distribution 
leads to enhanced counterion localization near the macromolecular surface. 

Oosawa-Manning (OM) condensation [2, 3] predicts that some of the counterions of a 
highly charged polyelectrolyte (PE) will "condense" near the surface of the PE, renormalizing 
the charge of the PE and leading to an effective Manning parameter £ e ff = 1. A rigorous 
definition of OM condensation can be derived from the PB solution for a homogeneously 
charged cylinder by taking the infinite dilution limit, p — > oo [21]. Above the "Manning 
threshold" £ = 1, a certain fraction / = 1 — l/£ of the counterions do not escape to infinity, 
and are therefore considered to be condensed; below the Manning threshold, there are no 
condensed counterions. For the two-state model, this definition of condensation needs to 
be modified slightly. In particular, all of the bound counterions in the TSM are clearly 
condensed. The free counterions, on the other hand, are still treated with the PB solution for 
a homogeneously charged cylinder; therefore, the results of OM theory still apply to them. 
In particular, if £ r > £oj then a fraction / = 1 — l/£ r of the free counterions are "Manning 
condensed" in the TSM, in addition to the bound counterions. Thus, when £ r > £0, the 
total number of condensed counterions in the TSM and in OM theory are in fact equal. It 
is important to distinguish between the bound and Manning condensed counterions here: the 
former are confined in a small volume near the surface of the PE, whereas the latter are 
contained within a more diffuse "Manning layer" in the neighborhood of the PE surface [21]. 
As is clear from eq. (4), it is also possible to have £ r < £o- In this case, no additional free 
counterions are Manning condensed in the TSM, and £ e // = £ r < 1. This can occur both when 
£ > 1, where there are some condensed counterions in OM theory, and when £ < 1, where OM 
theory predicts no condensation. In both cases, the fraction of condensed counterions in the 
TSM exceeds the prediction of OM theory. In other words, it is possible to have counterion 
condensation beyond that predicted by OM theory in the two-state model. It is important to 
note, however, that this is only possible at finite dilution; in the infinite dilution limit, it is 
straightforward to show that eq. (4) implies that the total number of condensed counterions 
in the TSM always agrees with OM theory. 

Comparison of the TSM with numerics. To test the validity of the TSM, we have 
solved the Poisson-Boltzmann equation numerically for a linear array of fixed charges. The 
discrete PB equation in this geometry is given by 




2^ B Q 2 r 2 (l- 7r cot[7 r ln(£)]) 2 



(5) 



PB[y(x)\ = 



1 



V 2 y(x) + 



N s e~ y{S) e{r - 5) 1 



]T S rfi S z , Zi - (6) 



4ttI b Q 



J <Pxe-yW6{r - S) 7r(Ar/2) 2 Az 



=1 



M.L. HENLE et al: COUNTERIONS AT DISCRETELY CHARGED SURFACES 



289 




2468 10 2468 10 



rlh r/l B 



Fig. 1 - a) Schematic representation of a linear array of charges surrounded by neutralizing counteri- 
ons. b-c) Fraction of enclosed counterions (f>(r) for 8 — 0.44, R = 104, and various values of Q and d. 
Symbols correspond to the full PB solutions, solid lines to predictions of the TSM, and dotted lines 
to the PB solutions for the homogeneously charged cylinder, b) Results for Q — 3, d = 6h (triangles) 
and Q = 3, d — 24 (squares), corresponding to £ = 0.5 (below Manning threshold) and £ = 1.5 (above 
Manning threshold), respectively. The fraction of bound counterions in the TSM is = 0.08 and 
0.13, respectively, c) Results for Q — 4, d = 84 (circles), corresponding to £ = 0.5. Here, /(, = 0.52. 



where 9(x) = 1 for x > and for x < , and zi are the locations of the surface charges, which 
are assumed to be monovalent for simplicity (i.e. Q s = 1). Due the azimuthal symmetry of the 
problem, the Laplacian V 2 = <9 2 + d r /r + <9 2 . The second derivatives are approximated using 
a symmetric three step finite difference equation. An explicit, forward time relaxation scheme 
is used to find the solution to eq. (6): y^ t+1 \x) — y^(x) — AtPB[y^(x)]. Convergence takes 
approximately 10 5 timesteps with At = 0.1. Physical lengths are determined by setting the 
lattice spacings Ar = Az = 0.11b- The boundary condition at r — R is dy/dr = 0. Periodic 
boundary conditions are used in the z direction, with enough lattice points in each period 
to contain 5 fixed charges. In fig. 1 we compare the function <f>(r) calculated using the TSM 
with the numerical solutions to the PB equation. To illustrate the effects of surface charge 
discretization, we also show the PB solution for the corresponding homogeneously charged 
cylinder. Vb, which does not need to be specified in the numerics, is set to 47r<5 3 /3 for the 
TSM. When the fraction of counterions in the bound state /& = 1 — £ r /£ is small, then the 
agreement between the TSM and the numerics is excellent (fig. lb). When is large, however, 
the TSM predicts a much stronger localization of the counterions than is seen in the numerical 
solutions (fig. lc). This is due to the fact that the PB equation does not include the discrete 
nature of the counterions; rather, the counterion distribution is a smooth function. That is, 
the possibility of surface charge-counterion pair association - an effect explicitly included in 
the two-state model - is excluded from the numerical solutions to the PB equation. When the 
binding energy of a surface charge-counterion pair is small - in other words, when Q is small or 
5 is large - the number of surface charge-counterion pairs is small, and the PB equation should 
accurately capture the effects of surface charge discretization on the counterion distribution. 
Thus, our numerical results provide verification of the two-state model in this limit. When 
the number of surface charge-counterion pairs is large, however, the PB equation is deficient, 
and we do not expect good agreement between the TSM and our numerical results. 

In conclusion, we have introduced a simple two-state model to describe the counterion 
distribution for discretely charged macromolecules. However, both Poisson-Boltzmann theory 
and strong-coupling (SC) theory [5,22] can be applied to the same systems as the TSM, and it 
is important to understand in which regimes each theory is valid. PB theory is valid when the 
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interaction of discrete surface charges with discrete counterions is sufficiently weak. SC theory, 
on the other hand, is valid in precisely the opposite limit, where the interactions between 
counterions can be treated as a perturbation to the surface charge-counterion interactions. 
Finally, in the TSM, interactions between counterions are fully included (albeit at the mean- 
field level), and the possibility of strong countcrion-surfacc charge interactions is accounted 
for by allowing discrete counterions to bind to the surface. Thus, we expect that the TSM 
is valid in regions between the SC regime and the PB regime, where the interactions of the 
counterions with both themselves and the discrete surface charges are important. This can 
be tested by performing Monte Carlo simulations, which we plan to do in a future work. 
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